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ABSTRACT 

The  requirements  for  failure  prediction  in  ceramics  using  ultra- 
sonics have  been  examined.  These  show  that  the  absolute  predic- 
tion of  failure  at  acceptable  stress  levels  and  lifetimes  requires 
high  frequencies,  in  the  50-400  MHz  range.  The  ability  to  utilize 
such  high  frequencies  for  flaw  detection  studies  is  primarily 
dependent  on  the  attenuation  of  the  material  in  this  frequency 
range.  Attenuation  studies  performed  on  a variety  of  ceramic 
materials  have  shown  that  the  coarse  grained  or  porous  ceramics 
are  excessively  attenuating,  whereas  the  fine  grained  fully  dense 
ceramics  are  acceptable. 

An  analysis  of  attenuation  using  numerical  scattering  cross  sections 
and  microstructural  parameters  has  demonstrated  that  the  attenuation 
is  entirely  predictable  from  the  large  extreme  of  the  micro structure. 
The  analysis  has  also  suggested  that  attenuation  measurements  may 
permit  the  Implementation  of  statistical  failure  prediction  at  high 
levels  of  confidence.  In  materials  that  are  not  amenable  to  absolute 
failure  prediction  using  high  frequency  ultrasonics. 
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1.0  INTRODUCTION 

A primary  restraint  on  the  structural  application  of  high  strength  ceramics 

Is  the  large  variability  In  their  strength.  The  distribution  In  fracture  strengths 

Is  related  to  the  size  distribution  of  defects  Introduced  during  the  fabrication 

and/or  finalizing  operations.1  A more  extensive  and  reliable  utilization  of 

ceramics  could  be  achieved  If  the  components  containing  the  most  deleterious 

flaws  could  be  Identified  and  hence,  eliminated  as  a source  of  premature  failure. 

The  most  deleterious  flaws  in  ceramics  are  not  uniquely  related  to  their 
1 2 

size.  * The  flaws  In  ceramics  at  the  large  size  extreme  are  usually  pores,  in- 
clusions or  large  grains  (the  latter  are  sites  for  machining  flaws  when  they  are 
in  close  proximity  to  the  surface).  Pores  and  inclusions  are  the  major  source 
of  premature  failure  in  fine  grained  (<  10  um)  ceramics,  while  large  grains  are 
an  equally  important  origin  of  failure  in  the  coarser  grained  [7  30  um)  materials. 
However,  the  role  of  an  inclusion  as  an  Initiator  of  failure  depends  on  Its 
nature,  as  well  as  its  size.  Some  are  highly  deleterious,  such  as  an  inclusion 
with  a small  elastic  modulus  (E)  and  thermal  expansion  coefficient  (o  ) compared 
to  the  matrix  or  an  inclusion  that  can  react  chemically  with  its  environs  to 
cause  a degradation  of  the  surrounding  matrix.  Others  are  relatively  innocuous, 

such  as  Inclusions  with  a similar  E,  *Q  to  the  matrix.  Failure  prediction  tech- 
niques that  rely  on  the  direct  determination  of  the  size  of  defects  must  therefore 

be  cognizant  of  the  failure  heirarchy  as  it  pertains  to  the  defect  and  material  of 

concern . 

It  is  the  initial  objective  of  the  present  study  to  determine  the  material 
systems  for  which  the  ultrasonic  technique  is  capable  of  obtaining  sufficient 
information  about  the  defect  size  and  composition,  that  the  absolute  prediction 
of  structural  failure  can  be  achieved  with  high  confidence.  Thereafter,  alternate 
failure  prediction  techniques  will  be  explored  for  those  materials  systems  which  are 
not  amenable  to  high  frequency  ultrasonic  analysis. 
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2.0  ULTRASONIC  REQUIREMENTS 


The  fracture  of  ceramics  usually  occurs  by  the  direct  extension  of  pre- 
existing flaws.  The  largest  sharp  crack,  am>  that  can  be  tolerated  in  a com- 
ponent that  is  expected  to  survive  for  a time,  t^,  after  inspection  is  given 

2 f 1 2/ (n-2 ) 


by1; 


m 


■c  2 At1aaZVZ(n-2)  J 


0) 


where  Kc  is  the  stress  intensity  factor  for  mechanical  extension  of  the  flaw, 
is  the  tensile  stress  (often  the  principal  tensile  stress  in  the  element  con- 

a 

taining  the  flaw),  Y is  a parameter  that  depends  on  the  flaw  shape,  and  A and  n 
are  parameters  that  define  the  susceptibility  of  the  material  to  slow  crack  growth. 

Approximate  evaluation  of  am  for  four  structural  ceramics  (two  hot  pressed 
silicon  nitrides  and  two  hot  pressed  silicon  carbides)  obtained,  assuming  that 
10,000  hours  is  an  acceptable  inspection  interval  and  300  MPa  is  an  acceptable 
stress-sustaining  capability,  has  shown  that  all  defects  which  produce  flaws  with 
radii  larger  than  1 0-1 0d  vm  must  be  detected.  For  most  defects,  equating  the 
defect  size  to  the  flaw  size  affords  a conservative  prediction  of  failure;1  and 
hence,  a defect  detection  capability  based  on  the  flaw  sizes  derived  from  eqn  (1) 
would  be  quite  adequate.  The  principal  exceptions  are  inclusions  that  either 
have  a smaller  E,  aQ  than  the  matrix,  or  have  degraded  the  matrix  by  chemical 
reaction.1  The  detection  requirement  for  these  precursors  can  be  substantially 
more  stringent  than  suggested  by  eqn  (1);  and  a capability  for  detecting  flaws 
even  smaller  than  ~ 10  pm  may  be  required.  Further,  to  avoid  the  rejection  of 
components  with  relatively  innocuous  inclusions  and  to  ensure  the  absence  of 
highly  debilitating  inclusions,  an  effective  ultrasonic  technique  must  be  able 
to  establish  the  type  of  the  defect  being  probed,  at  least  within  a limited  range 
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of  possibilities. 

Present  ultrasonic  techniques,  which  only  use  the  reflected  intensity  as  the 

defect  indicator,  are  not  capable  of  providing  unambiguous  information  about  the 

3 4 

flaw  size  and  character  and  hence,  do  not  satisfy  the  above  requirements.  * At 
best,  a probabilistic  estimate  of  failure  is  provided  (Fig.  1),  such  that  some 
deficient  parts  are  accepted  (type  I errors)  and  some  satisfactory  parts  are  re- 
jected (type  II  errors  and  'false'  signals).  The  final  failure  probability  is 
prescribed  to  be  the  most  acceptable  value  consistent  with  the  economic  constraints 

4 

that  result  from  excessive  rejection  due  to  type  II  errors. 

This  situation  can  be  greatly  improved  by  adopting  more  pertinent  measures 
of  the  defect  scattering  characteristics.  The  amplitude,  AR,  of  an  ultrasonic 
signal  received  by  a transducer  after  it  has  been  scattered  by  a defect  of  radius, 

5 

r,  is  related  to  the  transmitted  amplitude,  Ay,  of  a plane  wave  by  ; 

G^n]  [exp[-2R  a(f)  ]]  (2) 

e,  v is  the  wave  velocity  in 
the  material,  R is  the  distance  between  the  defect  and  the  transducers,  S is  the 
scattering  at  the  defect,  ct  is  the  attenuation  of  the  material,  and  G is  the 
transducer  gain  (including  bond  loss,  ?(/)).  The  first  term  in  eqn  (2)  is  simply 
the  loss  caused  by  beam  spreading.  (This  loss  can,  of  course,  be  extensively 
reduced  by  the  use  of  a focussed  beam,  but  the  determination  of  the  relation 
between  AR  and  Ay  for  non-planar  situations  is  not  a well-posed  problem.) 

The  second  term  derives  from  the  defect  scattering.  The  scattering  of  plane 
waves  by  defects  has  been  the  subject  of  several  recent  studies, and  in  con- 
sequence, numerical  scattering  calculations  can  now  be  performed  for  any  ellipsoidal 
defect.  For  example,  the  frequency  dependencies  of  the  scattering  of  longitudinal 
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A j [/*9»r]  = f (4it)3/2  r2  ] [S(/‘r’9)]  [/GT(/) 


where  f is  the  frequency,  e is  the  scattering  angl 
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plane  waves  have  been  calculated  for  two  typical  defect  types,  spherical  voids 
and  spherical  WC  inclusions  In  Si^N^,  at  three  defect  radii  of  potential  interest 
(Fig.  2a, b).  The  two  principal  features  to  note  are  the  strong  dependence  of  S 
on  frequency  (~/  ) when  kr<l  (k  is  the  wave  number,  equal  to  2-ir/X ) , and  the 
marked  difference  in  the  detected  frequency  characteristics  for  the  void  and  the 
inclusion  at  kr>l.  These  features  indicate  that  the  maximum  scattered  amplitudes, 
as  well  as  much  important  Information  about  the  character  of  the  defect,  can  be 
obtained  by  operating  in  the  range,  1/3  7 kr  < 3.  For  flaws  in  the  size  range 
10-100  urn,  the  frequency  regime  of  greatest  interest  is  thus  10  to  400  MHz. 

Th®  third  term  is  the  transducer  gain,  which  includes  all  losses  associated 
with  t transfer  of  acoustic  energy  to  electrical  energy,  including  bond  losses 
at  the  interface.  These  losses  can  be  minimized  by  optimizing  the  design  of  the 
transducer,  the  bonding  system,  etc. 

The  fourth,  and  final  term  is  the  attenuation  of  the  material.  The  atten- 

O 

uation  normally  increases  with  increase  in  frequency  and  microstructural  size  , 
and  could  thus  be  an  important  constraint  on  the  utility  of  ultrasonic  flaw 
characterization  in  the  appropriate  high  frequency  regime.  The  measurement  and 
analysis  of  attenuation  is  thus  the  key  initial  requirement  for  determining  the 
viability  of  ultrasonics  as  a failure  prediction  technique. 
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3.0  A> iENUATION 

3.1  Measurement  Techniques 

Attenuation  measurements  require  that  allowances  be  made  for  the 
nature  of  the  radiation  pattern  of  the  sound  beam,  the  effect  of  the  geometry 
of  the  sample,  the  effect  of  the  boundary  conditions  of  the  sample,  and  the 
energy  absorbed  by  the  transducer  from  each  succeeding  echo.  These  difficulties 
were  overcome  at  low  and  high  frequencies  by  different,  but  closely  related  pro- 
cedures. At  low  frequencies,  the  following  procedure  was  adopted.  (1)  A buffer 
was  used  between  the  transducers  and  the  specimen  in  the  form  of  a water  bath 
(1-25  MHz),  a fused  silica  rod  (25-120  MHz),  or  a sapphire  rod  (>120  MHz).  This 
procedure  removed  the  transducer  from  the  data  analysis  and  thus  eliminated  the 
need  to  account  for  any  energy  loss  or  echo  distortion  in  the  transducer.  (2) 

An  on-line  computer  (Data  General  Eclipse)  was  used  in  conjunction  with  a fast 
gate  and  an  analog-to-digital  converter  (Biomation)  to  obtain  the  Fourier  trans- 
forms of  each  received  echo.  This  procedure  was  necessary,  especially  for  the 
broad  band  transducers,  since  the  attenuation  was  found  to  vary  rapidly  with 
frequency  and  thus  required  a separate  determination  at  each  (monochromatic) 
frequency  within  the  band.  (3)  The  amplitude  of  the  same  frequency  component  in 
each  echo  was  corrected  for  the  diffraction  incurred  during  the  echo's  travel 
path  to  eliminate  the  losses  due  to  beam  spreading.  (4)  The  effect  of  the  boundary 
conditions  of  the  sample,  i.e.,  the  impedance  mismatch  between  buffer  and  samole, 
was  taken  into  account  by  comparing  the  attenuation  for  the  case  of  the  buffer  on 
on  both  sides  of  the  sample  with  that  of  the  buffer  on  one  side  only.  This  pro- 
cedure allows  the  calculation  of  the  reflection  coefficient  which,  for  these 
materials  with  low  attenuation  values,  represents  a sizeable  correction.  These 
procedures  plus  some  care  in  sample  preparation  (parallelism  of  opposite  faces, 
and  a sample/buffer  diameter  several  times  that  of  the  radiation  beam)  were 
found  essential  for  meaningful  measurements.  Although  the  measurements  were 
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tedious  and  time  consuming,  at  least  the  corrections  (3)  and  (4)  could  be  incor- 


porated into  a computer  program  which  was  designed  to  give  the  value  for  the 


attenuation  at  each  frequency  for  just  a few  initial  data  inputs. 


The  approach  employed  at  high  frequencies  was  to  work  with  a pulsed 


rf  signal  (~0.1  s)  which  was  long  enough  that  it  could  be  regarded  a single 


frequency,  but  short  enough  that  reflections  from  the  front  and  back  surfaces 


of  the  sample  could  be  easily  distinguished.  The  amplitude  of  the  reflected 


signal  from  the  end  of  the  buffer  rod,  with  no  sample  in  place,  was  then 


measured.  Assuming  the  end  of  the  buffer  to  be  a perfect  reflector,  this  gave 


a reference  amplitude.  A,  for  the  incident  wave  on  the  end  of  the  rod,  from  which 


the  effects  of  transducer  loss,  and  diffraction  and  attenuation  in  the  rod,  were 


removed.  An  optically  polished  sample  with  two  parallel  faces  was  then  pushed 


against  the  buffer  rod  with  a gold  foil  placed  between  the  rod  and  the  sample. 


The  amplitude  of  the  reflected  wave  from  the  interface  was  now  measured  to  be 


TA.  Hence,  the  amplitude  B of  the  wave  entering  the  sample  was 


B = A(1  - Irl2)1* 


Assuming  an  attenuation  o per  unit  length  in  the  sample  of  length  i,  the  signal 


returning  to  the  interface  would  be  B exp(-2at).  Hence,  the  resultant  echo 


reaching  the  transducer  should  have  an  amplitude 


C = B(1  - M2)**  exp(-2afc) 


It  follows  that 


= 1 ^n(l  - | r | 2 ) + tn  £ 


Since  the  first  return  gives  an  independent  measure  of  r,  it  is  possible  to 


measure  a.  Typically,  when  a sapphire  rod  is  used  to  contact  the  ceramic. 


there  is  a fairly  good  impedance  match,  so  most  of  the  incident  power  passes 
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into  the  ceramic  and  an  accurate  measurement  is  not  difficult  to  make.  An  error 

2 

analysis  can  be  performed.  We  find  that  for  errors  <S(|r|  ),  6 (A/C) , the  error  in 
a,  5(o),  is; 


*<5(a) 


lilii!)  + 

1 - 1 rl 2 


6(A/C) 

A/C 


(6) 


Thus,  the  technique  should  only  give  serious  errors  when  |r|  -*-1;  i.e., 
when  the  reflection  at  the  interface  is  very  large.  In  this  case  the  system  is 
better  employed  in  a more  coi.ventional  manner  to  compare  the  first  echo  from  the 
back  surface  of  the  sample  with  a triple  transit  echo. 

3.2  Test  Data 

Attenuation  data  have  been  obtained  for  MgO,  PZT  and  ZnS,  Si ^4  and 
SiC  (Figs.  3,  4 and  5).  The  MgO  was  a fully  dense,  hot 

pressed  material  with  an  equiaxed  grain  structure  (Fig.  6a)  and  no  noticeable 
porosity.  The  ZnS  was  also  fully  dense,  but  had  a more  complex  microstructure 
(Fig.  6b)  that  contained  many  twins.  The  silicon  nitride  materials  were  two 
fully  dense  hot-pressed  samples  (prepared  from  a-Si^N^  powder  with  2 and  5 percent 
MgO  additives)  and  a reaction  bonded  material.  The  coarsest  features  in  the  hot- 
pressed  materials  (Fig.  7a)  were  equiaxed  particles  of  silicides,  carbides,  etc. 
up  to  10  urn  in  diameter,  with  the  remainder  consisting  of  fine  elongated  grains 
of  silicon  nitride  up  to  ~ 5 pm  in  length  (Fig.  7b);  while  the  primary  micro- 
structural  features  in  the  reaction  bonded  materials  were  pores  (Fig.  8a, c)  and 
inclusions  (Fig.  8b)  up  to  60  pm  in  diameter.  The  silicon  carbide  samples  were 
made  from  a hot  pressed  material  and  a sintered  material.  The  hot  pressed 
material  contained  a proportion  of  relatively  large  equiaxed  inclusions  (Fig.  9a) 
and  a fine  (~  4 pm)  equiaxed  grain  structure  (Fig.  9d ) ; whereas  the  sintered 
material  contained  some  fine  porosity  (Fig.  9b),  and  tabular  grains  (Fig.  9c). 
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3.3  Analysis 

The  scattering  of  ultrasonic  waves  by  grains  has  been  presumed  by  pre- 
8 9 10 

vious  analysts  ’ ’ to  derive  from  the  difference  in  the  elastic  properties  of 
the  grain  and  matrix  (small  density  differences  have  been  neglected).  Then,  for 
a single  large  grain  in  a homogeneous  matrix,  the  Lame's  constants,  y and  X,  be- 
come the  important  scattering  parameters;  specifically,  the  scattering  cross 
section,  n,  is  determined  by  (aC)  /Cq,  where  C is  equal  to  X+2y  and  AC  is  the 
difference  between  the  modulus  of  the  grain  and  that  of  the  homogeneous  matrix, 
CQ.  However,  since  the  attenuation  is  the  sum  of  the  scattering  cross  sections 
of  a large  number  of  grains,  the  fundamental  scattering  cross  section  of  a grain 

in  a random  polycrystal  can  be  characterized  (for  attenuation  analysis)  by  the 

2 * 2 
average  elastic  mismatch,  < (aC)  >.  Further,  for  cubic  materials,  < (aC)  > can 


be  explicitly  written  in  terms  of  the  elastic  stiffness  constants,  C,  as; 


.8,10 


<(AC)2>  = (Cn  - C]2  - 2C44)2/75 


(7) 


8 10 

Equivalent  expressions  are  available  for  other  crystal  symmetries. 

The  numerical  value  of  the  scattering  cross  section  can  be  derived 

from  the  elastic  mismatch  by  determining  the  power  extracted  from  the  incident 

wave  by  the  scattered  dilational  and  shear  waves.  In  the  Rayleigh  regime  (kr<l), 

8 10 

analytic  expressions  have  been  derived,  and  for  an  incident  dilational  wave;’ 

5i 

(8) 


n = 


, 3,4,.  2 , 

4,1  < V-  (AC)2 


5 v, 


T 


(kr  « 1 ) 


= K V, 


In  textured  materials,  an  average  weighted  by  the  predominant  orientation  must 
be  used. 
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Note  the  strong  dependencies  of  the  scattering  on  the  grain  volume,  Vg,  the 
frequency,  /,  and  the  longitudinal  wave  speed,  v^.  At  shorter  wavelengths, 
kr  ***  1,  the  scattering  cross  section  must  be  determined  numerically.  This  can 
be  achieved  if  the  scatterer  is  described  in  terms  of  an  average  elastic  stiff- 
ness constants,  and  C^,  and  the  matrix  in  terms  of  the  polycrystalline 

elastic  properties.^  Then  the  energy  scattered  as  dilational  and  shear  waves 

12 

can  be  computed  using  a recently  developed  exact  scattering  theory.  The  fre- 
quency dependence  of  the  scattering  cross  section  determined  in  this  fashion  for 

large  grains  in  MgO  is  plotted  in  Fig.  10,  using  the  normalized  cross  section, 

2 

fl/2irr  , and  kr  as  coordinates. 

The  scattering  by  pores  can  be  derived  in  a similar  manner,  except  that 

the  elastic  moduli  of  the  scatterers  are  now  zero.  The  analytic  expression  in 

12 

the  Rayleigh  regime  is 


t 


n 


3 4 2 

4*  rvD 

— ~T~  gc 


(9) 


where  Vp  is  the  pore  volume  and  gc  is  a parameter  that  depends  on  the  elastic 
properties  of  the  host,  and  typically  ranges  from  5 to  30.  Computed  values  of 
the  scattering  cross  section  for  pores  in  Si^N^  are  plotted  in  Fig  10. 

It  is  apparent  from  the  strong  size  (V)  dependence  of  the  scattering 
cross  section  (eqns.  8,9)  that  the  scattered  energy  derives  primarily  from  the 
large  size  extreme  of  the  distribution  of  scatterers.  Hence,  if  the  largest  grain 
(or  pore)  in  a small  volume  element,  aV,  is  regarded  as  the  scatterer,  and  the 
surrounding  grains  within  that  element  are  considered  to  be  the  homogeneous 
matrix,  the  total  scattered  energy,  Uy,  (neglecting  interactions  between  adjacent 
scatterers)  is  the  sum  of  the  scattering  cross  sections  of  the  extreme  value 
grains  (or  pores)  over  all  volume  elements.^  Hence,  the  attenuation  (for  a given 
material)  becomes, 


10 
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i = 1/AV 

a (/)  = ^-uT(/)  = \ ni(v1,/)  (10) 

i = 1 
or, 

00 

a (/)  * \ / n(V.)fi(Vi,/)dV1 
o 

where  0^  is  the  scattering  cross  section  of  the  ith  scatterer,  aV  is  the  median 

AV  in  unit  volume  of  material,  is  the  volume  of  the  i^  scatterer  and  n(V^) 

is  the  number  of  scatterers  per  unit  volume  in  the  size  range  Vi  to  + d . 

The  problem  of  selecting  the  largest  value  from  a fixed  sample  size,  N,  is  an 

1 3 

extreme  value  problem.  Hence,  the  size  distribution  of  the  largest  scatterers 

should  be  characterized  by  one  of  only  three  possible  extreme  value  distribution 

functions  (Appendix).  This  substantial  restriction  of  the  possible  distribution 

function  greatly  simplifies  the  attenuation  analysis.  The  most  convenient  sample 

size  would  be  that  size  below  which  a consideration  of  additional  scatterers  pro- 

★ 

duces  no  significant  increase  in  attenuation. 

The  distribution  of  extreme  value  scatterers,  derived  from  chord  length 
measurements  (Appendix),  is  conveniently  represented  for  present  purposes  as  a 
histogram,  which  indicates  the  number  of  scatterers  per  unit  volume  in  a size 
(diameter)  range  D to  D+aD.  Grain  size  data  for  MgO  are  plotted  in  Fig.  11; 
while  pore  size  data  for  reaction  bonded  Si are  plotted  in  Fig.  12. 


★ 

The  effect  of  increasing  the  sample  size  in  a population  of  a fixed  total  size 
is  simply  to  distort  the  distribution  at  small  values  of  the  variate;  the  dis- 
tribution at  large  values  of  the  variate  is  essentially  unchanged.  Therefore, 
properties  determined  by  the  large  value  extreme  are  not  strongly  dependent  on 
the  sample  size.  The  minimum  sample  size  should  be  that  required  for  each 
large  value  scatterer  to  have  a complete  nearest  neighbor  array  of  matrix 
grains,  i.e.  3 in  a linear  element  or  ~ 9 in  a volume  element.  The  optimum 
sample  size  for  present  purposes  is  in  the  range,  ~ 5 to  10  for  linear  samples 
or  ~ 30-200  for  volume  samples. 


wr u Wf 
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The  attenuation  is  now  obtained  by  taking  the  scattering  cross  section 
from  Fig.  10  for  the  median  grain  (or  pore)  diameter  (ff)  In  each  range  D to  D+aD, 
and  summing  the  scattered  energy  over  the  total  number  of  scatterers  per  unit 
volume  (Figs.  11,  12)  within  this  size  range.  The  attenuations  predicted  for  MgO 
and  reaction  bonded  silicon  nitride,  plotted  In  Figs.  13a,b,  are  In  good  accord 
with  the  attenuation  data  for  both  materials.  Plots  of  the  contribution  to  the 
total  attenuation  from  each  size  range  of  grains  and  pores  (Fig.  14a, b)  confirm 
that  the  attenuation  derives  primarily  from  the  large  size  extreme,  and  shows 
that  the  predicted  values  should  be  slightly  underestimated  (as  observed,  Fig.  13) 
because  there  is  a small  contribution  to  the  attenuation  from  grains  (pores)  out- 
side the  measured  range  (especially  at  the  large  extreme).  Also,  as  might  be 
anticipated  from  the  scattering  cross  section,  the  attenuation  maximum  displaces 
to  smaller  sizes  as  the  frequency  increases  (Fig.  14). 


3.4  Implications 

The  success  of  extreme  value  statistics  and  scattering  cross  section 
computations  for  attenuation  prediction,  has  important  Implications  for  the 
analysis  of  attenuation  in  complex  (multi -component)  microstructures;  a problem 
which  has  previously  been  intractable.  The  large  extreme  of  each  component  within 
a microstructure  can  be  easily  evaluated  by  superimposing  a grid  on  a micrograph 
and  obtaining  extreme  value  chord  data  for  the  pertinent  phase  (Appendix).  The 
contributions  to  the  total  attenuation  from  each  microstructural  constituent  can 
then  be  assessed  from  the  scattering  cross  section.  The  dominant  attenuating 
phase  should  then  be  apparent.  This  information  would  suggest  microstructural 
design  approaches  for  minimising  attenuation,  which,  if  feasible,  would  enhance 
the  utility  of  high  frequency  ultrasonics  as  a failure  prediction  technique. 


12 
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4.0  DISCUSSIONS 

4.1  Flaw  Detectability 

The  relative  detectability  of  defects  in  various  ceramic  materials 
can  be  deduced  from  eqn  (2)  by  combining  the  calculated  values  of  the  intensities 
scattered  from  defects  and  the  bond  scattering  with  the  measured  material 
attenuation.  Typical  amplitudes  at  the  receiver  (relative  to  the  transmitted 
amplitude)  are  plotted  in  Fig.  15  for  spherical  voids  located  1 cm  from  the 
surface  in  hot  pressed  and  reaction  bonded  silicon  nitride.  For  the  hot  pressed 
material  the  peak  signal  amplitudes  are  obtained  above  50  MHz,  and  substantial 
amplitudes  are  retained  to  much  higher  frequencies;  whereas,  for  the  highly 
attenuating  reaction  bonded  material  the  peak  signal  amplitudes  are  depressed 
to  frequencies  below  20  MHz,  and  very  low  signal  amplitudes  obtained  above  40 
MHz.  High  frequency  ultrasonic  analysis  could  thus  be  viable  for  the  fully 
dense,  fine  grained  materials  that  exhibit  minimal  attenuation;  but  is  not 
expected  to  be  of  any  value  for  the  highly  attenuating  porous  or  coarse  grained 
materials. 

4.2  Role  of  Attenuation  Measurements  in  Failure  Prediction 

Although  ultrasonics  cannot  be  used  to  identify  discrete  defects  in 
materials  with  a coarse  microstructure,  it  is  suggested  herein  that  the  atten- 
uation provides  information  about  the  large  extreme  of  the  microstructure  that 
might  allow  the  statistical  prediction  of  failure  to  be  performed  with  acceptable 
confidence. 

Statistical  methods  have  previously  been  unacceptable  because  they 
make  the  unwarranted  assumption  that  the  flaw  population  from  batch  to  batch  is 
invariant.^  However,  ultrasonic  attenuation  measurements  on  specific  components, 
can  provide  information  about  the  flaw  population  in  that  component,  whenever 
the  flaws  are  inherent  flaws  that  relate  to  extreme  value  features  of  the 
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microstructure.  This  information  can  then  be  utilized  to  specify  the  failure 
probability  for  that  component  with  high  confidence. 

The  attenuation  at  a specific  frequency  is  related  to  the  distribution 
of  grain  (or  pore)  diameters,  D$,  by  (c.f.  eqn  10) 

oo 

2a  = j n(D$)fl  (D$)dDs  (11) 

o 

Hence,  by  applying  reasonable  functional  forms  from  n(D)dD  (such  as  the  second 

extreme  value  function  described  in  the  Appendix)  and  combining  with  the  computed 

n(0$)  for  the  scatterer  of  concern,  self  consistent  relations  between  a and 

n(Ds)dD$  can  be  obtained  by  performing  a few  attenuation  measurements  over 

the  pertinent  range  of  microstructural  variability.  The  probability  of  fracture 

1 14 

P(S)  of  a small  volume  element  AV..  at  an  applied  stress  level  S is  given  by;  * 


r 


l 


P(s)  = AVi  j n(Df ) d Df 
DfC($) 


(12) 


where  Df  are  the  diameters  of  the  fracture  controlling  flaws  and  Dfc  is  the  critical 
flaw  diameter,  at  the  stress  S, given  by; 


The  probability  of  fracture  D(Sm)  of  a component  of  volume  V subjected  to  a 

variable  stress  S(x,y,z), which  has  a maximum  value  Sm,  can  be  derived  from  the 

14 

product  of  the  elemental  fracture  probabilities  as; 

' i = V/AV.  • 0 

P(Sffl)  » 1 - exp  - ^ AVi  j n(Df  )d  Df  (14) 

. 1 = 1 Df  c ( S ) 
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Ir*  t 

*. 


Hence,  If  the  microstructural  elements  that  cause  scattering  are  the  same  as 
those  that  control  fracture,  n(D^)dD^  In  eqn  (14)  can  be  replaced  by  n (D$ )d  D$, 
derived  as  above  from  the  attenuation  measurements,  and  then  eqn  (14)  can  be 
solved  uniquely  for  the  fracture  probability,  P(Sm) , of  the  component.  There 
are  many  examples  of  coarse  grained  or  porous  ceramics  wherein  fracture  is 

15 

dominated  by  the  large  grains  or  pores,  e.g.  reaction  bonded  silicon  nitride, 
coarse  grained  oxides. This  approach  might  thus  be  tenable  for  materials  of 
this  type. 

Finally,  we  note  that  if  the  flaw  size  distribution  is  defined  by  the 
second  extreme  value  function  (eqn  A5),  then  at  the  very  large  extreme  that 
causes  fracture  it  reduces  to 

n(Df)  = (15) 

t 

which  on  substitution  into  eqns  (14)  and  (13)  gives  the  failure  probability; 

(16) 

where  SQ  is  a parameter  containing  k,  6,  Kc  and  V.  This  probability  distribution 
is  commonly  found  in  strength  tests  (the  ' two -parameter  Wiebull'  distribution), 
indicating  that  the  exponent  k that  enters  the  attenuation  behavior  is  related 
to  the  'Wiebull  modulus'  m,  that  determines  the  strength  variability,  by; 

m = k+4  (17) 

This  interesting  inference  merits  further  investigation. 


P(Sm)  - 1 - exp 


i - V/AV 


' v/a¥  ^ k+4 

S,  »-  / (if  « 
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5.0  APPENDIX 

EXTREME  VALUE  GRAIN  SIZE  DISTRIBUTIONS 

Grain  diameter,  D,  distributions  can  be  obtained  directly  from  planar  sections. 
The  simplest  approach  is  to  measure  the  lengths  of  chgrds,  L,  that  correspond  to 
grain  boundary  intersections  at  linear  traces  across  the  plane,  and  to 
apply  the  relation 

where  n(L)  is  the  measured  number  of  chords  per  unit  length  in  the  size  range  L 

to  L+dL  and  n ( D)  is  the  calculated  number  of  grains  per  unit  volume  in  the  size 

range  D to  D+dD.  Information  about  distributions  in  specific  grain  diameter 

regimes  can  also  be  estimated  in  this  fashion,  because  the  characteristics  of 

grainswithin  a specified  size  range  derive  from  the  corresponding  characteristics 

18 

of  the  chords  in  the  neighborhood  of  that  size.  Additionally,  for  extreme  value 

analysis,  since  the  fraction,  £,  of  a linear  trace  occupied  by  extreme  value  chords 

18 

is  also  the  fraction  of  the  volume  occupied  by  extreme  value  grains,  a fixed  sample 
size  used  for  chord  length  determinations  should  also  correspond  to  an  approximately 
fixed  number  of  grains  per  sample.  Hence,  the  extreme  value  nature  of  the  sampling 
is  essentially  retained,  and  the  extreme  value  chord  distributions  should  afford 
a close  (although  not  exact)  representation  of  the  characteristics  of  the  extreme 
value  grains.  The  extreme  value  distribution  can  thus  be  estimated  by  measuring 
the  lengths  of  the  largest  chords  in  equal  size  samples  along  a series  of  linear 
traces.  Then,  the  number  of  chords  per  unit  length  of  the  trace,  within  a 
specified  size  range,  yields  n(L).  Thereafter,  n(D)  can  be  obtained  from  n(L) 
using  eqn  (Al).  This  latter  procedure  is  facilitated  if  the  extreme  value  function 


describing  n(L)  can  be  obtained  from  the  chord  length  data,  such  that  eqn  (Al ) can  be 


t 
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transformed  by  performing  the  differential.  This  Is  achieved  by  plotting  the 
chord  length  data  in  terms  of  the  cumulative  probability,  «(L),  that  the  chord 
length  is  larger  than  L ($  is  determined  by  simple  order  statistics,  * = i/n+1), 
as  shown  in  Fig.  16  for  the  data  obtained  on  MgO.  Then,  the  functional  form  of 
<5 (L)  can  be  deduced  from  the  best  fit  extreme  value  data  plot.  For  most  single 
phase  materials,  it  has  been  found  (Fig.  16)  that  the  second  extreme  value  function1 


pertains,  namely; 


$(L)  = 1 - exp  - (Lq/L) 


where  Lq  and  k are  distribution  parameters.  Differentiating  eqn  (A2)  to  obtain 
the  frequency  distribution,  $(L),  gives 


♦ (L)  = 


But  n(L)  is  related  to  <J>(L)  by; 


-(Lo/L)k 


n(l)  ■ ^ 


where  L is  the  average  length  of  the  extreme  value  chords.  Substituting  n(L) 


from  eqn  (A4)  into  eqn  (Al)  and  performing  the  differential  thus  gives 

2k  £ L exp[-(L  /L)k]  \ / L \k  1 6 exp[-(L  /L)k] 

n(D)  = 2 [ 2+k-k  \ ~r~i  = 

* r(l-l/k)  Lk  ^ 1 v L ' J Lk+j 


n(D)  = 2+k-k  (^)  ] 

ir  L 


The  number  of  grains  per  unit  volume  with  diameters  between  L and  L+dL  can  thus 

be  obtained  directly  from  n (L ) and  the  parameters  k and  LQ  (the  latter  being 

derived  from  a best  fit  of  the  cumulative  probability,  *(L)  to  the  second  extreme 

we 

value  distribution).  F1nallyy(note  that  by  changing  variables  in  eqn  (A5),  n(2a) 
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reduces  at  large  values  (L  > LQ)  to  another  extreme  value  function  of  the  second 
type  (c.f.  eqns  (A3),  (A4) ) ; 

B “ 

n(°)  = y “5+T  exPt-(B0/e)°]  (A6) 

6 

3 

where  a = k/3,  8 * L and  and  y are  constants  containing  L , k and  £. 

This  result  confirms  the  earlier  statement  that  the  extreme  value  sampling 
of  the  chords  is  equivalent  (except  at  small  L)  to  an  extreme  value  sampling  of 
the  grains.  This  result  is  an  essential  part  of  the  extreme  value  approach. 


In  fact,  because  k is  large  (>5)  for  most  grain  and  pore  size  distributions 
(see  Fig.  16),  L does  not  need  to  be  much  larger  than  L for  eqn  (A6)  to 
pertain.  0 
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7.0  FIGURE  CAPTIONS 


Fig.  1.  The  relation  between  flaw  size  and  detection  probability  based  on 

signal  amplitude  for  two  ultrasonic  systems  or  sensitivity  settings. 

Fig.  2.  Calculated  scattering  from  defects  in  silicon  nitride,  (a)  spherical 
cavities,  (b)  spherical  WC  inclusions. 

Fig.  3.  Attenuation  data  for  three  relatively  coarse  grained,  fully  dense 
ceramics,  MgO,  ZnS,  PZT. 

Fig.  4.  Attenuation  data  for  two  hot  pressed  and  a reaction  bonded  silicon 
nitride. 

Fig.  5.  Attenuation  data  for  a hot  pressed  and  sintered  silicon  carbide. 

Fig.  6.  Reflected  light  optical  micrographs  of  etched  surfaces  of  (a)  MgO, 

(b)  ZnS. 

Fig.  7.  Micrographs  of  the  hot  pressed  Si-jN^  with  5%  MgO,  (a)  reflected  light 
micrograph  showing  inclusions  of  silicodes,  carbides,  etc.,  (b)  a 
scanning  electron  micrograph  of  an  HF  etched  fracture  surface  showing 


p 


the  elongated  grain  structure. 

Fig.  8.  Scanning  electron  micrographs  of  fracture  surfaces  of  the  reaction 
bonded  silicon  nitride. 

Fig.  9.  Micrographs  of  the  silicon  carbide  materials,  (a)  an  optical  micrograph 
of  the  hot  pressed  material,  (b)  an  optical  micrograph  of  the  sintered 
material,  (c)  a polarised  light  optical  micrograph  of  the  sintered 
material,  and  (d)  a scanning  electron  micrograph  of  a thermally  etched 
hot  pressed  material. 

Fig.  10.  Scattering  cross  sections  for  pores  in  silicon  nitride  and  grains  in 

magnesium  oxide;  r is  the  radius  of  the  scatterer  and  k is  the  wave  number. 
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Fig.  11.  A histogram  of  chord  lengths  (L)  and  grain  diameters  (D)  for  the  MgO 
shown  In  Fig.  6(a). 

Fig.  12.  A histogram  of  chord  lengths  (L)  and  pore  diameters  (D)  for  the  reaction 
bonded  Si-jN^  shown  In  Fig.  8. 

Fig.  13.  Predicted  attenuations  compared  with  the  measured  values  for  (a)  MgO 
(b)  reaction  bonded  SijN^. 

Fig.  14.  The  predicted  dependence  of  the  attenuation  on  the  size  of  the  scatterer 
at  several  frequencies,  (a)  MgO,  (b)  reaction  bonded  Si^. 

Fig.  15.  Predicted  relative  detected  amplitudes  for  hot  pressed  silicon  carbide 
and  reaction  bonded  silicon  nitride. 

Fig.  16.  Chord  length  data  for  MgO  fitted  to  the  second  extreme  value  distribution 
function. 
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